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ABSTRACT 

We present Keck/HIRES data with model atmosphere analysis of the 
helium-dominated polluted white dwarf GD 40, in which we measure atmo- 
spheric abundances relative to helium of 9 elements: H, O, Mg, Si, Ca, Ti, Cr, 
Mn, and Fe. Apart from hydrogen whose association with the other contami- 
nants is uncertain, this material most likely accreted from GD 40's circumstellar 
dust disk whose existence is demonstrated by excess infrared emission. The data 
are best explained by accretion of rocky planetary material, in which heavy el- 
ements are largely contained within oxides, derived from a tidally disrupted 
minor planet at least the mass of Juno, and probably as massive as Vesta. The 
relatively low hydrogen abundance sets an upper limit of 10% water by mass in 
the inferred parent body, and the relatively high abundances of refractory ele- 
ments, Ca and Ti, may indicate high-temperature processing. While the overall 
constitution of the parent body is similar to the bulk Earth being over 85% by 
mass composed of oxygen, magnesium, silicon and iron, we find n(Si)/n(Mg) = 
0.30 lb 0.11, significantly smaller than the ratio near unity for the bulk Earth, 
chondrites, the Sun, and nearby stars. This result suggests that differentiation 
occurred within the parent body. 

Subject headings: planetary systems - stars: abundances - white dwarfs 
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1. INTRODUCTION 

Detailed information about the bulk compositions of planetary bodies is generally lim- 
ited by the lack of available mechanisms to probe planet interiors. However, the technique 
of measuring accreted element abundances of tidally disrupted planetary bodies that pol- 
lute the atmospheres of cool white dwarfs provides an indirect, yet sensitive method to 
study bulk compositions. The otherwise pure white dwarf atmospheres are a clear canvas 
for displaying the relative bulk elemental abundances of the swallowed planetary parent 
body/bodies, as demonstrated with Kcck/HIRES observations, and subsequent abundance 
analyses, of the dramatically polluted white dwarf GD 362 (Zuckerman et al. 2007, Koester 
2009a, Juraet al. 2009a). 

In contrast to main sequence stars whose atmospheres contain mainly primordial el- 
ements, cool (6,000 K < Tgff < 25,000 K) white dwarfs provide a relatively clean slate 
since their high surface gravities cause heavy elements to settle downward, leaving above 
a nearly pure hydrogen and/or helium atmosphere (Schatzman 1945). The settling times 
for "high-Z" (atomic #Z > 2) elements to diffuse down from the outer convection zone are 
always short compared to the white dwarf lifetime (Paquette ct al. 1986), and thus the 
high-Z elements cannot be primordial. The competing outward force of radiative levitation 
is negligible below 25,000 K (Chayer et al. 1995), and convective carbon dredge-up, which 
may play a role in bringing carbon to the surface for Tg// < 13,000 K (Koester et al. 
1982), cannot explain the presence of heavier elements. Consequently, high-Z atmospheric 
contaminants in cool white dwarfs must be accreted from an external source. Given the 
purity of white dwarf (WD) atmospheres, even trace amounts of external pollutants can be 
detected spectroscopically, especially with high resolution, high signal-to-noise ratio (SNR) 
observations. 

While heavy element absorption features have been detected in spectra of single WDs 
for over half a century (Kuiper 1941; see also Sion et al. 1983; Wesemael et al. 1993 
for current WD classifications and references), it is relatively recent that the growing 
number of discoveries of infrared excess emission (due to light-reprocessing by the dust of a 
circumstellar disk) at many of the most heavily polluted WDs has elucidated the nature of 
the accretion reservoirs (Zuckerman &: Becklin 1987; Becklin et al. 2005; Kilic ct al. 2005, 
2006; von Hippel et al. 2007; Jura et al. 2007a, 2007b; Brinkworth et al. 2009; Farihi et al. 
2009; Melis et al. 2009, in preparation). The spectral energy distribution of the infrared 
excess is well modeled by a disk geometry and the dust temperatures typically range from 
inner to outer edges between 1200 K and 300 K (Jura et al. 2007a), which correspond to 
radii between 10 and 100 R* (0.13 - 1.5 R©), within the tidal radius of a typical white 
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dwarf (Jura 2003) . Double peaked Ca II emission line profiles at some of the hotter (T^f f 
near 20,000 K) WDs reveal the dynamics of gaseous material in Keplerian rotation also 
within 1.5 Kq (Gaensicke et al. 2006, 2007, 2008). Given the location, the dust and 
gas must have been produced (or delivered) after the stars' red giant phase. Additional 
insight into the nature of the circumstellar material comes from the lO/xm silicate emission 
features detected with Spitzer IRS spectra of dusty white dwarfs (Jura et al. 2007b, 
2009b, Reach et al. 2005, 2009), which can be modeled by micron-sized olivine-like glasses 
that produce a characteristic emission profile resembling that of circumstellar, but not 
interstellar grains. Evidence is strong that these disks are derived from the tidal disruption 
of an asteroid or minor planet (see Jura 2008) which has undergone orbital perturbation 
(Dcbcs & Sigurdsson 2002) during post-red giant evolution of the star's planetary system. 

The material in the atmospheres of cool white dwarfs with detected disks has almost 
certainly been accreted from their circumstellar reservoirs. The atmospheric contaminant 
abundances are linked to the composition of the accreted source material (Koester 2009a) . 
This is potentially a rich source of highly detailed information regarding the nature of 
extrasolar planetary matter. In their work on GD 362 with Keck/HIRES, Zuckerman et 
al. (2007) detected 15 high-Z elements, a major advance over the low resolution discovery 
work that detected Ca, Mg, and Fe (Gianninas et al. 2004). Zuckerman et al. (2007) found 
element abundances consistent with a refractory-rich, volatile-poor composition which they 
conclude very probably derived from an asteroid or minor planet. What is needed are 
more high-spectral-resolution, high-SNR observations and abundance analyses of additional 
WD/perturbed-planetary systems. In this paper we present a detailed analysis of the 
atmospheric composition of GD 40 (= WD 0300-013), an exemplar of polluted and dusty 
white dwarfs. For the first time, ground-based observations of an externally polluted WD 
have detected all the major constituents of terrestrial planets: O, Mg, Si, and Fe. It 
is a useful reminder that these 4 elements compose 94% of Earth's mass. With current 
models and calculations of the relationship between accreted and atmospheric abundances 
in polluted white dwarfs (Koester 2009a, Jura et al. 2009a), we focus on understanding 
the bulk chemical composition of GD 40's polluting parent body. We show that rocky 
extrasolar planetary bodies exist (ed) in this presumably ancient planetary system and 
may have experienced differentiation during formation and/or evolution. 



2. OBSERVATIONS & MEASUREMENTS 



The data were acquired in good weather conditions with the high resolution echelle 
spectrometer, HIRES (Vogt et al. 1994) on the Keck I telescope at Mauna Kea Observatory. 
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A 1.148 arcsecond slit width was used for all observations, providing a spectral resolution 
of R = A/AA ~ 40,000. Table 1 lists observation dates and setups. 

GD 40 (Vmag = 15.56, spectral type DBAZ) is the prototype of the DBZ class, being 
the first helium-hned cool white dwarf observed to display heavy element traces (Wickra- 
masinghe et al. 1975, Shipman &; Greenstein 1983). HST FOS observations in the UV 
(Shipman et al 1995), and the subsequent analysis by Friedrich et al. (1999) and Wolff et 
al. (2002), hereafter F99W02, yielded abundance measurements of H, C, Mg, Si, Ca, and 
Fe, with uncertainties ~ 0.5 dex. An optical spectrum of GD 40 (Koester et al. 2005a) 
was obtained in the ESO Supernova Progenitor (SPY) survey, and has been used in de- 
termining the effective atmospheric temperature, Tg// = 15,300 K, with log g set to 8.0 
(Voss et al. 2007), which we adopt for the models used in this work (§ 3 contains further 
discussion of T^ff and log g). Recently, GD 40 was found to possess an infrared excess 
(Jura et al. 2007a) and lO^m silicate emission of circumstcUar origin (Jura et al. 2009b). 
These discoveries made Keck/HIRES observations of GD 40 a priority. 



2.1. Data Reduction 

Data reduction was carried out with standard IRAF^ tasks for echelle data. A PyRAF 
code meta-structure controls the data processing flow, maintains file organization and adds 
versatility. The data have undergone typical echelle reduction procedures: bias removal, 
flat-fielding, variance weighted ("optimal") extraction, dispersion calibration with lineariza- 
tion and heliocentric doppler correction. Much of the echelle blaze function was removed 
by dividing each spectrum by a low-order polynomial fit to the extracted spectrum of the 
quartz lamp flat exposure. It was desirable to increase the signal-to-noise ratio by combin- 
ing orders of overlapping wavelength coverage, and to reduce the known HIRES artifacts 
(waves and bends in the continuum of a wavelength scale greater than ^10 A) most likely 
due to variable vignetting (Suzuki et al. 2003). Since we do not require absolute flux 
calibrated information, and we are mainly interested in heavy-element absorption lines 
(equivalent widths < 3 A), the echelle orders were continuum normalized to unity with a 
low-order polynomial (typically 5th order cubic spline) to reduce artifacts and align the 
echelle orders. An error associated with this action was estimated by varying the order 
of the normalization polynomial, which resulted in a < 1% variation in equivalent width 
measurements of typical heavy element absorption lines. Adjacent orders of overlapping 



^http://iraf.noao.edu/docs/recommend.html; Churchill (1995) is also instructive 
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wavelength coverage were combined in an average, weighted by the blaze function to sup- 
press the contribution from the noisier edges of the CCD chip (e.g. Aoki 2006^), and a 
final order-merged, normalized spectrum was obtained. 

We employed different instrumental setups to achieve our desired wavelength coverage 
and sensitivity. For the 2008 red data, an additional re-normalizing processing step had 
to be applied to calibrate and remove 2nd (diffraction) order blue flux contamination in 
the region A8200-9000A. The 2006 red data had a different instrumental setup, providing 
an un-contaminated comparison to check our renormalization procedure. The final all- 
combined CD 40 HIRES spectrum has SNRs in the range of 35 to 50 for the spectral 
interval 3300 - 4500A, and between 60 - 90 for 4500 - 7000A. Shorter than 3300A, the 
SNR is ~ 20; near 7000A, the SNR is ~ 50; redward of SOOOA, the SNR is < 35. 

2.2. Spectral Lines 

CD 40's optical spectrum is dominated by large (equivalent widths ~ 1-10 A) pres- 
sure broadened helium lines, a signature of the spectral class "DB". The high-resolution 
and high-SNR Keck/HIRES spectrum displays absorption lines of many more species. In 
addition to He we identify the presence of 9 elements, H, O, Mg, Si, Ca, Ti, Cr, Mn, and 
Fe (Figures 1-8), whose detections are securely established through multiple (> 4), un- 
ambiguous lines for all elements heavier than hydrogen. Hydrogen is detected via a weak 
but clearly visible Ha feature, also observed by Voss ct al. (2007). Except for the Ca II 
H&K fines and two Ti II fines (3242.0 A and 3383.8 A), CD 40's absorption lines arise 
from excitation levels above the ground state, and arc therefore not of interstellar origin. 
The very large and broad profiles of the Ca II H&K fines indicate their atmospheric origin. 
The carbon detection from UV data (Shipman et al. 1995, F99W02) brings the number of 
known elements aside from He in GD 40's atmosphere to 10. 

Each absorption line was visually inspected. Line identifications made use of the Ku- 
rucz database^, and the Vienna Atomic Line Database. We used IRAF (splot) to fit 
the line profiles in the continuum-normalized spectrum with the line center, core depth, 
FWHM, and the local continuum value as free fit parameters. The equivalent width (EW) 
and a fitting error were calculated in each profile fit. All of the lines analyzed were best 
fit with a Lorentzian profile and the fit-derived EWs were checked for agreement with 



^http://optik2. mtk.nao.ac.jp/~waoki/QL/specana200601e.pdf 
^http://www.pmp. uni-hannover.de/cgi-bin/ssi/tcst/kurucz/sekur. html 
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direct numerical integration. Three different measurements were performed for each line, 
varying the interval of nearby continuum data points from which the continuum level was 
determined. The reported EW is the average of the 3 measured EW values, and the total 
EW uncertainty includes the RMS of the 3 fit values added in quadrature with the average 
fitting error and the 1% error from continuum normalization (§ 2.1). Measurement uncer- 
tainties range from ~ 3% - 50%, with the peak of the distribution at ~ 10%, depending 
on the strength of the line and the ability to establish a well-defined continuum. Our line 
identification and measurement results are listed in Table 2. Prom 99 high-Z absorption 
lines, the average measured (heliocentric) radial velocity of GD 40, which includes its grav- 
itational redshift, is 10 km s~^, with a standard deviation of 2 km s~^ (Figure 9). The 
average high-Z line velocity agrees with the average helium line velocity (9 ± 4 km s~^) 
implying a coincident physical position with helium in the gravitational well, an additional 
confirmation of atmospheric location for the absorbers. 

Shipman et al. (1995) tentatively identified lines of neutral oxygen in the 1 A reso- 
lution UV HST (FOS) spectrum, noting a possible blend with silicon. However, Friedrich 
et al. (1999) re-analyzed those spectra to derive abundances, and attributed the uncertain 
features entirely to Si I and Si II. Our ground-based optical/near-IR detection of oxygen in 
GD 40 clears up the previous ambiguity regarding this element. Oxygen is an important 
element for study since it is the most abundant constituent of rocky planetary material 
and is a component of water. 

Nickle is not unambiguously detected in the HIRES data, but two of the stronger 
expected transition lines are possibly associated with noisy features (e.g. A3465.6 in Figure 
5). Aluminum was also tentatively identified previously in the FOS spectrum (Shipman 
et al. 1995), and fit for abundance by Friedrich et al. (1999), although the latter authors 
regard its detection as unconvincing. 

3. ELEMENT ABUNDANCES 

The basic method of determining elemental abundances in a stellar atmosphere is to 
compare the observed spectrum with a synthetic spectrum derived from a model of the 
star's atmosphere. The input physics of our models are described in Koester (2009b). Ab- 
sorption line data were taken from the Kurucz database, the Vienna Atomic Line Database, 
and Melendez & Barbuy (2009). Key input parameters are the effective temperature, T^ff-, 
surface gravity, g, and the atmospheric abundances of elements relative to the dominant 
constituent, in this case, helium. At GD 40's temperature, along with a very weak hydrogen 
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feature, the strong helium hnes indicate a hehum-mass-dominated atmosphere. 

Voss et al. (2007) derived the effective temperatures and surface gravities, of a large 
sample of helium-line-dominated (DB) WDs from the broadening of helium lines in the 
SPY spectra. For GD 40, they found T^ff = 15,300 ± 500 K with log g set to 8.0 (cooling 
age = 0.2 Gyr, mass = 0.59 Mq). In DB WDs with effective temperatures below ~ 16,500 

K, log g cannot be determined unambiguously from the helium line broadening, so Voss 
et al. (2007) adopted the typical value, log g = 8.0. Since 97% of DB WDs with effective 
temperatures between 16,500 and 19,000 K have log g in the range 7.7 to 8.3 (Voss et al. 
2007), we assume a ±0.3 dex uncertainty in log g. 

Our synthetic spectrum reproduces the observed absorption spectrum (see Figures 1 - 
8), and we compare the EWs of the model lines with those measured in the data. The EW 
measurements are translated into abundance values relative to helium by interpolation 
between models with different relative abundances. In a representative curve-of-growth 
plot of the Si II A3862.6 line in GD 40 (see Figure 10), the EW versus abundance relation 
is well defined. For each element, the adopted abundance value is the weighted average 
of the set of abundances from individual lines of that element. We choose to weight the 
contribution from each line by the strength of its detection, i.e. using weights = SNR at 
line center (column 4 of Table 2) = [line core depth] /[RMS nearby continuum]. This gives 
nearly the same result as weighting by the inverse square of propagated measurement er- 
rors. For each set of high-Z clement lines, the standard deviation of the mean (SDOM) of 
the set of abundances are in the range 0.03-0.09 dex, which always exceeds (by 50-200%) 
the measurement error propagated from EW uncertainties. This is expected since system- 
atic errors enter in going from EW to abundance. We take the SDOM as the "formal" 
error of the abundance average. To be conservative in estimating our errors we include the 
difference between weighted and unweighted averages as an additional contribution to the 
uncertainty. Finally, we varied the model parameters by their uncertainties, Tg// (15,300 
lb 500 K) and log g (8.0 it 0.3), one at a time, to calculate a per-element abundance error 
associated with the model temperature and gravity uncertainties, respectively. In Table 3 
we list the derived abundances and uncertainties, relative to helium, for the 9 detected con- 
taminants and we place useful upper limits on Al and Ni, as well as 6 additional elements. 
For the abundances relative to helium, the T^ff uncertainty is usually the dominant source 
of error, as can be readily seen from Table 3. 

With a precision of ~ 0.1 - 0.2 dex, our work reduces the abundance uncertainties 
for previously identified elements, Mg, Si, and Ca, by at least a factor of 5. For Mg, Si, 
Ca and Fe, our new abundances are all lower than the nominal values in F99W02 by ~ 
30-40 %, which for Mg, Si, and Ca are within the errors (see Figure 11). Differences in 
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model temperatures used for this analysis (15,300 K) and F99W02 (15,150 K) are not 
the cause of the systematic abundance offset, since at lower Tg// the implied abundances 
are even smaller. The HIRES detection of Ha (EW = 0.23 A) agrees with the SPY data. 
However, based on an Hq EW measurement of 1.0 A (Greenstein & Liebert, 1990), F99W02 
obtained a hydrogen abundance nearly an order of magnitude higher than the HIRES and 
SPY results. In retrospect, the feature assigned to Ha (see Figure lb of Greenstein & 
Liebert, 1990) was noise. Prom their possible UV detection of Al, Priedrich et al. (1999) 
suggest log[n(Al)/n(He)] = -7.15, just below our detection limit of log[n(Al)/n(He)] < -7.1. 
In our composition analysis we include the small, but relevant contributions from Ni and 
Al using the HIRES-derived upper limits on their abundances. 

In principle, Teff can also be evaluated from the abundance agreement between differ- 
ent ionization states of an clement, since we expect a consistent value for the total number 
of atoms of the element. Both Mg and Pe display lines from two different ionic states, 
Mg I/Mg II and Pe I/Pe II. The separated abundances for the different ionic states and 
varied model input temperatures are shown in Figures 12 and 13. Although the mean 
deviations of the abundances are close to the differences between the average values, it 
is possible to distinguish between models with different temperatures. For both Mg and 
Pe, the ionization balance is improved at the lower temperature limit, and degraded at 
the higher temperature limit. Also, the ionized species always reflects a lower abundance 
than the neutral atoms, likewise suggesting the model temperature is too high. With our 
presentation of the uncertainties in Table 3, it is straightforward to recover all the element 
abundances, to within ~0.05 dex, that result from the model with Tg// = 14,800 K (if ever 
in the future this becomes established), from the Table 3 data by subtracting the Teff un- 
certainty from the abundances listed in column 2. Similarly, the 15,800 K abundances can 
be re-calculated by raising all the column 2 abundances up by the T^ff uncertainty. Such 
potential shifts in model temperature do not significantly alter our results of the relative 
elemental abundances in the parent body causing GD 40's atmospheric pollution. This is 
because the abundances of polluting elements all move together, either down or up, with 
varying Teff, the pollutant-to-pollutant ratios are less sensitive to the model temperature 
than are the pollutant-to-helium ratios. 



4. RESULTS AND DISCUSSION 



How do these abundances help us understand extrasolar planetary systems? 
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4.1. Bulk Composition 

First we note that GD 40 has accreted at least the mass of a large asteroid. With 
a helium-dominated convective zone (CVZ) that contains a fraction q = log {Mcvz/M) = 
-5.86 (Koester 2009a) of the white dwarf's total mass of 0.59 Mq, there is 1.7 x lO^^ g of 
helium in the CVZ. Prom this and Table 3 abundances, we find a total mass of pollutant 
material presently detected in the CVZ of 3.6 x 10^^ g (Table 4). For comparison, the 
mass of asteroid Juno is 3.0 x 10^^ g with a radius of 117 km (Konopliv et al. 2006), 
although Konopliv et al. (2006) note that Juno has a near metallic density, 4.43 it 0.43 g 
cm^'^ A 3.6 X 10^^ g asteroid of density 3.0 g cm~^ has a radius of 142 km. 

Second, wc identify the material to be in a category of terrestrial-planetary compo- 
sition. Figure 14 depicts the bulk composition of GD 40's pollutants and a variety of 
solar system objects for comparison. It is immediately apparent how GD 40's volatile-poor 
abundances are very different from primitive solar system abundances, and we find the 
closest agreement with an Earth-like composition. 94% of Earth's mass comes from just 
4 elements: O, Mg, Si, and Fe (Allegre et al. 1995). Ca, Ni and Al contribute another ~ 
5%, and all other elements are minor constituents. Similarly, by bulk mass, the material 
polluting GD 40's atmosphere is 88% O, Mg, Si, and Fe. Ca contributes another 6%. Up- 
per limits for the CVZ mass percentages of Al (2.5%) and Ni (1.3%) have been included 
in the totals. 

Jura (2006) and Farihi et al. (2009) have pointed out that, for a handful of polluted 
WDs with carbon measurements, including GD 40, the polluting abundance of carbon 
relative to iron is down by a factor of 6 or more from the solar ratio. The deficiencies of 

carbon and hydrogen arc consistent with rocky matter. By similar arguments as apply for 
GD 362, which have been described in detail by Zuckerman et al. (2007), based on the 
quantity and constitution of pollutants, an interstellar or cometary origin for the high-Z 
elements is "exceedingly unlikely". 

The origin of the hydrogen is less certain because unlike the high-Z elements, it does not 
settle below the outer convective zone, which results in an accumulation of all the primordial 
and accreted hydrogen over the WD lifetime. Both Voss et al. (2007) and Dufour et al. 
(2007) find a trend of increasing hydrogen mass with increasing age in helium-dominated 
WDs which argues strongly in favor of an accreted, rather than primordial, origin. GD 
40's hydrogen abundance is consistent with these populations of helium-dominated WDs 
(Figure 1 of Jura et al. 2009a). Since we do not have a constraint on its arrival epoch, 
it is possible that some, or maybe all, of the presently observed hydrogen is a result of 
previous accretion of other planetary matter (whose high-Z components have settled), or 
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due to accretion from the interstellar medium, and may not even be associated with the 
most recently accreted material. Alternatively, some, or possibly all, of GD 40's hydrogen 
could be tied to the currently observed heavy element pollution. At best we can place 
some upper limits. One path by which hydrogen may survive processes that result in a 
refractory-rich environment is through H2O. Jura &: Xu (2009) propose that water from 
accreted minor planets may have a role in the hydrogen "pollutions" of helium-dominated 
WDs. Our measured hydrogen abundance (at its upper error limit) imposes an upper 
limit to the amount of water or ice which could have been accreted (in total for GD 40's 
entire WD cooling time) at 3.6 x 10^^ g. If all of this water was accreted along with the 
present high-Z material, then an upper limit to the amount of water in the bulk pollutant 
composition is 10%. Wc show below, in §4.3, that due to settling of the heavier elements, 
the upper limit is likely to be much less than this. In either case GD 40 has not accreted 
a very large amount of water, or hydrogen, along with its high-Z material. 

GD 40's atmospheric pollution is consistent with having accreted rocky material. How- 
ever, the atmospheric abundances, n(Z)/n(Hc), need not be the same as the accreted abun- 
dances, and going from the measured abundances to the accreted ones requires some insight 
and modeling. 

4.2. Accretion — Diffusion 

The relationship between observed and accreted abundances is governed by the inter- 
play of accretion and diffusion. Key inputs are the characteristics of the encounter with 
the reservoir, the rate that matter is falling into the star and the rates by which heavy 
elements are settling (diffusing) downward. Fundamentals of accretion - diffusion theory 
have been described in detail by Dupuis et al (1993). More recently, Koester k, Wilken 
(2006) and Koester (2009a) have quantified mixing zone mass fractions and settling times 
of commonly accreted elements for a variety of white dwarf atmospheres in the temperature 
range of interest. Table 4 lists settling times, tset{Z) (= diffusion time, T(el) in Koester 
2009a), calculated specifically for GD 40's CVZ. 

Following the models of Jura ct al (2009a) and Koester (2009a), accretion is occuring 
from the observed disk of material and may be represented by an exponentially decaying 
mass function with a characteristic lifetime, i^isfc (= ^acc in Koester 2009a). We are ob- 
serving the system at a time tot^, which we would like to constrain, with respect to tget 
and tdisk- While exact relations governing the abundance ratios can be found from equa- 
tion 8 of Koester (2009a) and equation 8 of Jura et al. (2009a), here we summarize the 
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characteristics of three hmiting cases: (1) In the early phase of accretion, toi,s ^ tgetj 

the 

accreted abundance ratios are directly the observed atmospheric ones. (2) For times on 
the order of the diffusion time scales a steady state is approached, and if the reservoir is 
long-lived, tget ^ tobs ^ tdiski the steady state values of the accreted abundances are 
modified from the measured abundances by the ratios of the settling times (equation 7 of 
Koester 2009a, equivalently equation 16 of Jura et al. 2009a): 



The ratios of settling times are at most a factor of 2.6 between detected elements 
studied in this work. (3) After the accretion reservoir is spent, tdisk < iobs, the atmospheric 
element ratios deviate exponentially with time from the accreted ratios, as gravitational 
settling causes the heavier elements to sink more rapidly than the lighter ones. This case 
leads to a significant departure from the accreted abundances. 

Referring to GD 40 CVZ abundances shown in Figure 15, the high-Z ratios do not 
differ from the Earth or CI chondrite ratios by more than a factor of 4, (except for carbon), 
mostly the factor is 1-2 as is the case for O, Mg, Ti, Cr, Mn and Fe. This is not at all the 
exponential deviations between observed and accreted abundances that are expected from 
gravitational settling if accretion has ended, paused, or significantly decreased. Also, there 
is no trend of increasing depletion of element abundances (relative to Earth or chondrite 
abundances) with decreasing element settling times. Thus, the relative atmospheric abun- 
dances are consistent with accretion being presently ongoing, which is further supported 
by the observational existence of a disk. This implies that we are either observing GD 40 
in a steady state, at an early phase of accretion, or somewhere in between. Figures 14 and 
15 depict these two limiting extremes of the accreted abundances. 



The Earth-like composition of GD 40's pollutants is almost certainly not coincidental, 
but rather, it is probably linked to the stoichiometric balance of rock-forming oxides. In 
the solar system rocky planetary bodies are primarily composed of mineral oxides (MgO, 
AI2O3, Si02, CaO, Ti02, Cr203, MnO, FcO, Fe203, NiO, and minor constituents) and 
possibly metals, e.g. Fe, Ni, Cr, etc, particularly if a core has formed. 




(1) 



4.3. 



Oxides 



It is a remarkable result that the quantities of heavy elements in GD 40's atmosphere 
are consistent with a rocky composition wherein heavy elements are largely contained 
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within oxides. In other words, the observed oxygen is of sufficient quantity to form minerals 
with the other observed elements, as is evident from computing the number of oxygen atoms 
needed to form oxides, 2^(^)0^(2), (measured relative to He): 

n{0) oxides _ q{Z) n{Z) 
n{He) ~ ^ p{Z) n{He) 

With Table 3 observed element (Z) abundances, which directly represent early-phase accre- 
tion as defined in §4.2, Equation 2 adds up roughly to 1.5 x 10~^, significantly below the 
measured oxygen abundance, n(0)/n(He) = 2.5 x 10^^. Additionally, as discussed in §4.1, 
some of the oxygen may be associated with water/ice, since a portion, or possibly all, of 
the atmospheric hydrogen could be coeval with the current accretion episode. Constrained 
by the relatively low hydrogen abundance, at most 0.5 x 10~^ of the observed n(0)/n(He) 
can originate from water. 

Now we have the question of accounting for all the excess oxygen, since the refractory- 
rich, volatile-poor element abundances imply that unbound, volatile oxygen is not a likely 
vehicle. What happens in the steady state? In order to calculate the steady state abun- 
dances from Equation 1, pollutant-to-pollutant ratios are needed. We list early phase and 
steady state abundance ratios, relative to oxygen, in Table 5. To evaluate the oxygen bud- 
get, let's look at the oxide balance in greater detail. As is suitable for a rocky parent body, 
we assume oxygen is mainly carried in oxides or water. Removing the helium dependence 
of Equation 2 and working with the oxygen ratios, a balanced oxygen budget requires: 

1 n{Hpar) sr^ q{Z) n{Z) 

2 n{0) ^p{Z)n{0) 

We use Rpar to represent the undetermined fraction of the observed hydrogen that is 
associated with water in the parent body causing the current accretion. Another ambiguity 
comes from the fact that there are three principle states of iron in the Earth: metallic Fe (as 
in Earth's core), FeO, and Fe203, but in the white dwarf's atmosphere we are measuring 
the sum of Fe atoms from all states. Since we do not have a clear determination of the 
partitioning of Fe in the GD 40 parent body, in Table 5 we show sums for 3 different 
distributions of Fe in order to better understand the range of effects on the oxide balance. 
Iron is potentially one of the main carriers of oxygen, and in order to account for the 
observed O abundance, we find that the Fe can not be all metallic. From the steady 
state oxide balance, we need at least 45% of the iron atoms in FeO (at least 30% if the 
more oxygen-rich species, Fe203, is the principle state) in order to carry the oxygen. This 



(2) 



(3) 
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constraint means that by total mass the accreted material was likely less than 28% metallic 
iron. 

Referring to the oxygen budget of Table 5, if one includes all of the hydrogen as water 
and assumes that all of the iron is in the more oxygenated form, Fe203 (no metallic Fe), 
then, even at the upper limit of the uncertainty, the fraction of oxygen accounted for by 
early accretion phase abundances gets close to, but does not quite add up to 1. In other 
words, the parent body abundances that would be associated with early phase accretion 
have too much oxygen compared to what can be carried by high-Z elements and water. 
A more likely scenario is that accretion has been going on long enough for the system to 
be in a steady state with heavier elements settling more rapidly than the oxygen. Indeed, 
we see in Tabic 5, the steady state abundances are in balance with respect to oxides (i.e. 
the oxygen budget adds up to unity), and in this configuration, the problem of too much 
oxygen in the parent body disappears. 

We learn from this exercise that the oxide balance favors the steady state accretion 
situation (represented by the steady state bulk composition in Figure 14 and the steady 
state abundances in Figure 15, with element-to-element ratios derivable from the steady 
state column of Table 5). In order for the system to achieve a steady state, the disk 
lifetime, must be longer than the settling time scales tdisk ^ ^set ~ 6 x 10^ yr. By these 
arguments, a lower bound on GD 40's disk lifetime is tdisk > 10^ years. With GD 40 
in, or near, an accretion - diff'usion steady state, then after just a few settling times, the 
total accreted mass would be a few x lO^^g, or about the mass of Vesta, the second 
most massive asteroid in our solar system (Konopliv et al. 2006). We have constrained 
tset < iobs < tdisk J but we do not know how many settling times have passed. Also, since 
the disk is opaque (Jura et al. 2007a) we cannot estimate its mass, so we do not know how 
much more mass may be waiting to accrete from there. 

Returning to the question of water content in the parent body, since we have argued 
in favor of accretion taking place for at least a couple of settling times, then in the steady 
state, an upper limit to the percent of water by mass is 3%, decreasing nearly linearly with 
tabs ^ tdisk (equation 6 of Jura et al. 2009a). The water content of a planetary body 
can be a distinctive measure. By mass, the Earth's mantle"^, Ceres and Callisto are 0.03 - 
0.1%, 25%, and 50% water, respectively (van Thienen et al. 2007, McCord & Sotin 2005, 
Canup &; Ward 2002, respectively). Meteorites originating from the outer asteroid belt 



*The light element composition of Earth's core is uncertain and may include an appreciable amount of 
hydrogen and oxygen. 
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(2.5-4 AU) are up to 10% water by mass, while those originating from the inner asteroid 
belt (~ 2 AU) range from 0.05-0.1% water by mass (Morbidelli et al. 2000). 

4.4. An Extrasolar Minor Planet 

GD 40 has a disk within its tidal radius, circumstellar silicate dust emission, the mass 
of a large asteroid in its convective zone, and pollutant abundances consistent with rocky 
material. The most plausible origin of the atmospheric heavy elements is a tidally disrupted 
planetary parent body, such as an asteroid/minor planet. Presumably, this object was part 
of an ancient planetary system, had its orbit perturbed (Debes & Sigurdsson 2002), such 
that it journeyed within the star's tidal radius and was disrupted (Jura 2003), forming a 
disk of material that is accreting onto the star. 

In Figure 15 the element abundances are shown normalized to the primitive material CI 
chondrites to assist our interpretation of how the GD 40 parent body may have evolved from 
a primitive constitution. We see an overall trend of increasing abundances with increasing 
condensation temperature indicative of thermal processing. The close agreement with bulk 
Earth values for O, Mg, Cr, Mn, and Fe anchors the Earth-like abundance nature of the 
polluting body, permitting us to interpret the deviations, most notably the depletion of 
silicon and the build up of refractory elements, in terms of evolutionary processes. 

A distinctive result we find is n(Si)/n(Mg) = 0.30 it 0.11 in significant contrast to 
the ratio near unity for the bulk Earth, the Sun, chondrites, and nearby stars. This is not 
part of the pattern of the thermal trend, nor can it be caused by differential gravitational 
settling, since Si and Mg have nearly the same diffusion times (see also Figure 15). Our 
result is consistent with the F99W02 value of n(Si)/n(Mg) = 0.3lo;2. In this paper, we 
have reduced the uncertainties over previous work enough to recognize that the primordial 
elemental mix in this evolved planetary system has been altered within the parent body. 

Figure 16 depicts the distinction from cosmic abundances. The GD 40 star and proto- 
planetary-system was probably not born with a peculiar Si/Mg ratio. Presently no solar 
system object has an observed Si/Mg abundance with such a low ratio. Perhaps the closest 
comparison to a measured material in our own solar system comes from the Pallasite mete- 
orites, which are thought to derive from the mantle-core boundary of a differentiated body 
(i.e. physically stratified into mantle, core, and/or crust), and measure n(Si)/n(Mg) = 0.6 
(Buseck 1977). One possibility is that the distinctive Si/Mg ratio may be associated with 
mantle/crust differentiation, in which Si was concentrated in a crust, and Mg in a mantle, 
and the crust was later lost. Evidence for eroding collisions exists in our solar system in 



-15- 



the form of the Moon and the diversity of the population of asteroids and meteorites such 
as the Pallasites. The so-called hit-and-run process (Asphaug et al. 2006), thought to play 
a role in terrestrial planet formation, results in increasingly exotic characteristics for the 
unaccreted remnants of planetesimal collisions. If differentiation takes place in a planetary 
body, Si-rich minerals tend to rise into the crust, whereas Mg-rich minerals prefer the man- 
tle. In Earth's crust n(Si)/n(Mg) = 11, but in a large rocky planetary body such as Earth, 
the 35 km crust is a relatively thin layer and the mantle/crust mass ratio is 169:1. Thus, 
even though Si is highly concentrated in the crust, a relatively small amount by mass has 
been depleted from the Earth's mantle. In an asteroid sized body, such as Vesta (radius 
= 270 km), the crust is thought to be ~27 km, and a mantle/crust mass ratio has been 
estimated at 3:1 (Ghosh & McSween 1998, Jones 1984), making a substantial Mg and Si 
differentiation in an asteroid plausible. 

A difficulty of this explanation is with the observed enhancements of Ca, Ti, and 
possibly Al in GD 40. Similar to silicon, the elements Ca, Ti and Al also tend to concentrate 
in the crust of a differentiated body. For instance. Si, Ca, Ti, and Al are respectively, 1.3, 
1.8, 3.4 and 3.5 times more abundant in Earth's crust than in Earth's mantle (Palme k, 
O'Neill 2007). If GD 40's parent body lost its crust, then shoTild not Ca, Ti and Al also 
be depleted along with Si? Something else must be going on. Since Ca, Ti and Al are 
refractory, most likely their enhancements arc attributable to a significant thermal process. 
Indeed, Earth's mantle, while depleted in Si relative to primitive CI chondrites, is modestly 
enhanced in Ca, Ti and Al (see Table 6), notably due to the refractory lithophile nature 
of these elements. Given the strong enhancements of Ca and Ti in GD 40's polluting 
planetary body, it likely experienced a more extreme thermal history than Earth, as may 
have come about for example from intense heating during the red giant phase or throTigh 
the energy of collisions. Another consideration is that exposure to high temperatures 
during planetary formation and/or evolution could be important. With an estimated WD 
mass of 0.59 M© for GD 40, its main sequence progenitor was likely greater than 1.5 M© 
(Weidemann 2000). As previously suggested by Jura et al. (2009a), the circumstellar 
environments around stars more massive than our Sun may support conditions of more 
intense heating of the protoplanetary material, planetesimals and planets than in the solar 
nebula, which could explain the dramatically high abundances of calcium polluting the 
white dwarfs GD 362 and GD 40. One possibility for GD 40 is that both heating and crust 
loss have occurred, and that Ca and Ti are enhanced because the heating effect on their 
abundances dominates over the crust loss. In that case, if the GD 40 parent body had not 
lost its crust, the Ca, Ti, and Al abundances would be even more enhanced than what are 
measured now. 
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Fegley & Cameron (1987) describe another mechanism - a silicate vaporization process 
- which removes substantial quantities of material (e.g. ~ 80% of Mercury's silicate) from 
the outer layers of a differentiated body due to prolonged heating. Their calculations were 
designed to explain the unusually high density of planet Mercury and are successful at 
doing so by modeling the protoplanet's exposure to a high-temperature (2500 K - 3500 K) 
phase of the primitive solar nebula. More generally we obtain a sense of the constitutional 
evolution of a planetary body with an initially chondritic composition after experiencing a 
period of intense heating. The models predict a unique resulting abundance pattern such 
that, relative to chondritic abundances, the mantle+crust becomes depleted in Si02, while 
MgO, CaO, Ti02, and AI2O3 are enhanced. In particular, model #4 (from Table 4 of 
Fegley &: Cameron 1987) yields Si/Mg, Ca/Mg, and Ti/Mg ratios in agreement with CD 
40's measured values, as shown in Table 6. We do not refer to the iron abundances from 
the mantle+crust results of Fegley & Cameron (1987) since in their model most of the 
iron resides in a core. This silicate vaporization model can simultaneously reproduce the 
Si depletion, and the Ca and Ti enhancements in CD 40's polluting parent body, however, 
the calculated high density planet has an iron-to-silicate mass ratio twice that of Earth, 
which is appropriate for planet Mercury, but is not what we observe in CD 40. Although 
this model involves differentiation of the parent body and essentially an erosion of the outer 
layers - mechanisms which we think are important - it is highly idealized and its application 
to real systems is uncertain. 

Finally we note that a crust/mantle differentiated body is likely to have also formed 
a core composed primarily of Fe metal as in Earth's core. As discussed in §4.3, within the 
uncertainties of CD 40's steady state pollutant abundances, the oxide balance allows up to 
70% of the iron to originate from a metallic form (< 55% if the iron oxide is predominantly 
FeO), although a rocky composition with no metallic iron is also allowed by the oxygen 
budget (Table 5). This translates to an upper limit on the size of a metallic core in the 
accreted parent body at ~25% of its total mass (assuming a Vesta-like core composition of 
0.89 Fe by mass, as modeled by Ghosh &; McSween 1998). For comparison. Earth's core is 
32% of its mass (Allegre et al. 1995), and estimates of core sizes in Mars, Vesta, and the 
Moon are 30%, 10% and 5% of their total masses, respectively (Righter k, Drake, 1996). 



5. CONCLUSIONS 

Our Keck/HIRES observations and abundance analysis make CD 40 the second dusty 
polluted white dwarf (CD 362 is the first) with a comprehensive set of measured element 
abundances and at least the mass of a large asteroid in its outer convective layer. CD 
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40's atmosphere is polluted by material with a rocky chemical composition, approximately 
similar to bulk Earth which is dominated by O, Mg, Si, and Fe. The inferred parent body 

was not a predominantly icy object, containing at most 10%, but more likely < 3%, water 
by mass. The balance of elements within oxides and water implies that the system is in, 
or near, an accretion - diffusion steady state, which places a lower bound on GD 40's disk 
lifetime at tdisk > 1 Myr. These results, taken together with the infrared Spitzer observa- 
tions, can best be explained by the accretion of an orbitally-perturbed, tidally-disrupted, 
ancient, planetary parent body of at least the mass of Vesta (radius ~ 270 km), carrying 
with it a history of this extrasolar system's formation and/or evolution. The suite of ele- 
ment abundances reveals enhanced Ca and Ti, probably due to thermal processing, and a 
distinctively low Si/Mg ratio which may be an indication that differentiation has occurred 
in this planetary system, since a minor planet stripped of its crust could conceivably cause 
the silicon deficiency. Within the uncertainties of the oxygen budget, the accreted body 
may have had a metallic core up to one-quarter of its total mass. Because a solid planetary 
body has conveniently been ground down to its elemental constituents, this dusty polluted 
white dwarf provides us a powerful tool for measuring the bulk chemical composition of an 
extrasolar minor planet, at a level of detail and precision greater than what is possible for 
most terrestrial planetary bodies in our own solar system. 
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Table 1: Observations 



UT Date 


Collimator 


A range 


exposure 






(A) 


(seconds) 


2006 Sep 09 


Red 


5,690 - 10,200 


2400 


2007 Nov 20 


Blue 


3,110 - 5,990 


2 X 2700 


2008 Nov 14 


Red 


4,500 - 9,000 


1800 


2008 Nov 15 


Red 


4,500 - 9,000 


2 X 1800 



Notes. The continuum SNR is discussed in §2.1 
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Table 2:: GD 40 Absorption Lines in HIRES Spectra 
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Table 2 — Continued 
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Cr II 


3118.6 


11.8 


4 


29 




5 


Cr II 


3120.4 


11.9 


4 


70 




30 


Cr II 


3125.0 


9.5 


7 


83 




7 


Cr II 


3128.7 


13.7 


2 


20 




5 


Cr II 


3132.1 


8.5 


8 


92 




6 


Cr II 


3147.2 


12.3 


3 


26 




14 


Cr II 


3180.7 


10.0 


5 


30 




10 
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Table 2 — Continued 
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11.0 


5 


37 




4 


Fe II 


3210.4 


10.6 


6 


54 




5 


Fe II 


3213.3 


11.6 


13 


105 




6 


Fe II 


3227.7 


10.8 


14 


141 




4 


Fe II 


3247.2 


9.5 


3 


33 




11 


Fe II 


3255.9 


6.7 


3 


28 




5 


Fe II 


3258.8 


8.5 


3 


32 




7 


Fc II 


3259.1 


10.7 


4 


28 




4 


Fc II 


3493.5 


11.2 


3 


17 
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Fe II 


4233.2 


8.2 


2 


50 




10 
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Table 2 — Continued 



Ion 






SNR^ 








(A) 


(km s"-*-) 


(line 

center) 


(mA) 




Fe II 


5018.4 


9.9 


7 


50 ± 


4 


Fe II 


5169.0 


9.2 


15 


81 ± 


3 


Fe II 


5276.0 


9.6 


2 


24 ± 


9 



Notes. - Prominent He I lines are detected at wavelengths (in A): 3188, 3820, 3889, 3965, 
4388, 4472, 4713, 4922, 5016, 5876, 6678 and 7281 (4026 and 7065 lie in gaps of the HIRES 
spectrum and were not accessible). The average He line velocity is 9 km s~^ with a stan- 
dard deviation of 4 km s~^. 
" Laboratory wavelength in air. 

^ Velocities are in a heliocentric rest frame; Ha has an asymmetric profile and slightly 
ambiguous line center (see also Figure 7); Mg II 4481 is a blended triplet - the line centers 
have not been individually determined here and the EW is the total for the blend. 

Line depth (from profile fit) divided by the RMS of the nearby continuum. 

§2.2 describes the equivalent width measurements and uncertainties. 
^ The 3154.2 line is an unresolved blend of Ti II + Fe II (Figure 2), not used in the abun- 
dance determination. 
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Table 3: GD 40 Atmospheric Element Abundances 
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Li 


iog[n(^Zj;/n(^±iejj 
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rormai weighted 
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11 


Si 


-6.76 -F0.07/-0.08 


7 1 


e 
o 


13 


Ca 


-6.88 -h0.14/-0.22 


D 4 


'i'i 


20 


Ti 


-8.61 -F0.13/-0.20 


7 1 


oO 


6 


Cr 


-8.31 +0.12/-0.16 


Q 1 


zy 


6 


Mn 


-8.62 -K0.13/-0.18 




ol 


4 


Fe 


-6.48 -F0.09/-0.12 


7 9 


zz 


1 




previous UV 








c 


-7.2 -I-0.2/-0.2 










1 1 T ^1 u >t' 1 1 in 1 1 ^ 


from ihc; a})S(;nce 


of (A/A) 




N 


< -5.3 


8680.3, 8216.3 






Na 


< -7.0 


5890.0, 5896.0 






Al 


< -7.1 


3961.5, 3944.0 






P 


< -8.2 


3509.0, 3786.6 






Sc 


< -9.8 


3572.5, 3613.8 






V 


< -9.0 


3118.4, 3125.3 






Ni 


< -7.7 


3407.3. 3465.6^, 3514.0^ 




Sr 


< -10.1 


4077.7, 4215.5 







Notes. - Abundances derived assuming a model with reyj=15,300K and log ^=8.0. 
The listed number abundance errors are the sum of contributing uncertainties added in 

quadrature. Carbon abundance from UV data (Friedrich et al. 1999, Wolff et al. 2002). 
"Formal" is the standard deviation of the mean for the set of element lines. For H with 
just one line the formal error is the measurement error. 

"Weighted" = half the difference between weighted and unweighted averages. 
"Tejj" and "log g" errors come from varying the model parameters by their uncertainties, 
re/y^=15,300ib500 K and log (/=8.0ib0.3. Since the element abundances all move either 
up or down together with varying temperature, the T^ff uncertainty directly represents 
this systematic shift, revealing the abundances that result from models with 14,800K and 
15,800K (see last paragraph of §3). 

''^ marginal features may be present at these wavelengths (e.g. Nill 3465.6 in Figure 5) 
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Table 4: Gravitational settling times, contaminant masses and steady state mass accretion 
rates for GD40's convective zone. 



atomic # 


Z 


log (iset/yrs) 


Mc,,/(102i g) 


M„ec/(10« g S-1) 


1 


H 


floats 


0.29 ± 0.11 




6 


C 


5.7943 


0.32 ± 0.19 


0.16 ± 0.10 


7 


N 


5.7720 


< 29.7 


< 15.9 


8 





5.7587 


16.6 ± 3.2 


9.2 ± 1.8 


11 


Na 


5.6680 


< 1.0 


< 0.7 


12 


Mg 


5.6623 


5.9 ± 1.9 


4.1 ± 1.3 


13 


Al 


5.6625 


< 0.91 


< 0.63 


14 


Si 


5.6520 


2.06 ± 0.35 


1.46 ± 0.25 


15 


P 


5.5986 


< 0.083 


< 0.066 


20 


Ca 


5.4694 


2.2 ± 0.9 


2.4 ± 0.8 


21 


Sc 


5.4389 


< 0.003 


< 0.004 


22 


Ti 


5.4373 


0.050 ± 0.018 


0.058 ± 0.021 


23 


V 


5.4398 


< 0.022 


< 0.03 


24 


Cr 


5.4321 


0.11 ± 0.03 


0.13 ± 0.04 


25 


Mn 


5.4105 


0.056 ± 0.019 


0.069 ± 0.023 


26 


Fe 


5.4006 


7.8 ± 1.9 


9.8 ± 2.3 


28 


Ni 


5.3715 


< 0.50 


< 0.67 


38 


Sr 


5.2089 


< 0.003 


< 0.006 


Total detected 




35.5 


27.4 



Notes. - Settling times, tget (= diffusion times, Tdiff) are defined in Koester (2009a). 
Contaminant masses are calculated from Table 3 abundances and a helium-dominated 

convection zone mass of Mc„2(He) = 1.7 x 10^'' g. 
In the steady state limit, Macc(Z) = Mcvzi'Z')/tset{'^)- 

Some upper limits, nitrogen in particular, are not very restrictive but are included 
nonetheless for completeness. 
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Table 5: Limits of Parent Body Abundances 



n 
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/r^(n\ 
n{L)/n{U) 


n{L)/n{U) 
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1 1 
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n n7i -I- n m fi 

U.U ( i ± U.UID 


n no -1- n no 
u.uy It U.U/ 
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U.UUiU It U.UUUZ 


n nn9i -i- n nnn/i 

U.UUZi It U.UUU4 


Pt 


n nnon ^ n nnn/i 

U.UUZU zL U.UUU4 


n nn/i9 -i- n nnns 

U.UU4Z zL U.UUUo 


Mn 


n nni n + n nnns 


n nn99 + n 0006 


Fe 


0.13 ± 0.02 


0.30 ± 0.05 


Ni 


< 0.008 


< 0.02 


Oxygen budget: 






H2O 


< 0.2 


< 0.1 


oxides (3 cases): 






all Fe in Fe203 


0.69 lb 0.09 


1.12 lb 0.12 


all Fe in FeO 


0.62 ± 0.09 


0.97 ± 0.10 


\ FcO. \ Fe luclal 


0.5G ± 0.09 


0.82 ± 0.09 



Notes. - Upper portion: two limiting extremes for the accreted abundance ratios relative 
to oxygen. Early phase and steady state are defined in §4.2. The pollutant-to-pollutant 
ratios are less sensitive to variations in Te// than are the poUutant-to- helium ratios, 
which results in smaller errors (30-60% smaller) listed here than what one obtains in 
propagated errors from Table 3. With the exception of hydrogen, the steady state values 
have been calculated using the settling times from table 4 according to Equation 1. For 
hydrogen, the steady state fraction relative to O depends on the number of settling times 
that have passed. An upper limit is estimated assuming tobs > 2 tset, i-e. n(H)/n(0) < 
(0. 28-1-0. ll)/2 (= early phase upper limit / 2 settling times). 

Lower portion: The oxygen budget, which includes Al and Ni at their upper limits, is 

the fraction of the observed oxygen that may be accounted for by the total sum of oxides 
(MgO, AI2O3, Si02, CaO, Ti02, Cr203, MnO, Fc in FC2O3, FeO or Fe metal, NiO) and 
water. Three different compositions involving iron are shown. Uncertainties from the 
element ratios are propagated in quadrature. The oxide, Na20 can also contribute, but we 
do not have an abundance measurement, or meaningful upper limit, for Na. In the silicate 
Earth, Na20 ~ Cr203 « 0.4% by mass, within the errors of our oxygen budget. The 
expectation is that the contribution from oxides and water will add up to unity (Equation 
3). The steady state abundances account for all of the observed oxygen, in agreement 
with a mineral stoichiometry. 
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Table 6: Distinctive Abundance Ratios 



Ratio 


Silicate 


GD 40 


GD 40 


F & C 


F & C 




Earth 


cvz 


Steady State 


model 4 


model 3 


Number ratios relative to CI chondrite values 


Si/Mg 


0.83 


0.31 ± 0.12 


0.32 ± 0.13 


0.28 


0.43 


Ca/Mg 


1.16 


4.0 ± 2.1 


6.2 ± 3.3 


3.0 


4.1 


Ti/Mg 


1.16 


1.8 ± 0.9 


3.1 ± 1.6 


2.6 


4.2 


Al/Mg 


1.16 


< 1.7 


< 1.7 


1.4 


2.0 



Notes. - According to Fegley &; Cameron (1987), mantle+crust vaporization of a 
planetary body due to exposure to high temperatures can result in the depletion of Si 
and the enhancement of Ca, Ti, and Al in the parent body. Fegley & Cameron data 
are from their Table 4. Their model 4 representing a nonideal magma heated to 3500 K 
predicts mantle+crust ratios in agreement with GD 40's pollutant abundances. Model 3 
for a nonideal magma heated to 3000 K is also near agreement with GD 40 abundances. 
Data for the silicate Earth (Earth's mantle+crust) are from McDonough k, Sun (1995); 
CI chondrites from Lodders (2003). 
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Fig. 1. — A portion of the Keck/HIRES spectrum of GD 40. Wavelengths are in air 
and a hehocentric rest frame. The model spectrum shown in red has been convolved 
with a Gaussian profile to reflect the instrumental broadening (resolution, R ~ 40,000), 
as measured by Gaussian fits to arc lamp lines. The model is velocity shifted from the 
laboratory rest frame by the mean value, 10 km s~^, measured from the data (see Figure 
9). The strongest Cr lines from Table 2 are displayed. 
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Fig. 2. — Similar to Figure 1, with prominent Ca lines that do not originate from the 
ground state. 
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Fig. 3. — Similar to Figure 1, with the strongest Fe hnes from Table 2 displayed 
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Fig. 4. — Similar to Figure 1, with the strongest Ti lines from Table 2 displayed 
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Fig. 5. — Similar to Figure 1, but with all Mn lines from Table 2 displayed. For ease of 
viewing all 4 Mn lines in one plot, these data are smoothed with a 7 point boxcar average, 
and the model resolution is matched to the smoothed data by convolution with a gaussian 
profile (which effectively corresponds to degrading the spectral resolution from R=40,000 
to R=20,000). 
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Fig. 7. — Similar to Figure 1, but for Ha. The line center in the data appears shifted from 
the systemic velocity (represented by the red model curve; discussed in §2.2) by ~ +15 km 
s-i. 
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Fig. 8. — Similar to Figure 1, but with all of the O lines from Table 2 displayed. The inset 
data of the 01 A8446 line are smoothed with an 11 point average. 
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Fig. 9. — Line velocities from Table 2. All heavy element spectral lines in GD 40 have 
unambiguous identifications and their measured line centers lie within a distribution of 
velocities centered on a mean heliocentric value of 10 km s~"^. 
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Fig. 10. — Example curve-of-growth for the Sill A3862.6A line in GD 40's atmosphere. 
The solid curve simply connects the points with a smoothed line. There is a high fidelity 
between the model equivalent width and input Si abundance. For a given absorption line 
the element abundance is determined by interpolation between models. 
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Fig. 11. — GD 40's atmospheric abundances from Table 3 derived from Keck/HIRES 
optical spectra compared with previous abundances derived from UV data (C, Mg, Al, Si, 
Fe) and optical data (both previous H measurements). Al and Ni are upper limits. 
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Fig. 12. — Magnesium abundances derived from different ionic states of magnesium witli 
models at low, central, and high values of GD 40's effective temperature. Error bars 
represent one standard deviation of the mean added in quadrature with the difference 
between weighted and unweighted averages for the set of lines of the given ion. The 
ionization balance is improved at the lower model temperature. 
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Fig. 13. — Same as Figure 12, but for iron. 
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100 




Fig. 14. — Bulk compositions. GD 40's volatile-deficient polluting abundances are most 
similar to a bulk terrestrial planet composition. Totals do not reach 100% because minor 
constituents are omitted. GD 40's mass total includes Al and Ni at their upper limits (not 
shown). Helium and sulfur are not measured in GD 40, but given their volatility, these 
elements are expected to be minor constituents. GD 40 data are from Table 5 for the two 
limiting extremes of parent body abundances, steady state and CVZ (= convective zone = 
early phase); Sun and CI chondrites from Lodders (2003); comet Halley from Lodders & 
Fegley (1998); bulk Earth from Allegre et al. (1995) 
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Fig. 15. — Abundances shown normalized to Mg and CI chondrite composition to empha- 
size the alteration relative to primitive abundances. Ni and Al are upper limits. O, Mg, 
Fe, Cr, and Mn are consistent with bulk Earth abundances. Similar to Earth, the volatiles 
C, O are depleted from primitive values. The refractory elements Ca and Ti are enhanced 
in the body that pollutes GD 40. The trend of increasing abundances with condensation 
temperature most likely is a result of thermal processing. The depletion of silicon is distinc- 
tive and suggests that some kind of differential processing is likely to have occurred in the 
parent body. The GD 40 data are from Table 5 (CVZ = convective zone = early-phase); 
Earth from Allegre et al. (1995) with Earth carbon data from Lodders & Fegley (1998); 
CI chondrites from Lodders (2003). 
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Fig. 16. — Data for 180 F and G dwarf stars from Reddy et al. (2003). For a sample of 
6 early B-type stars, Przybilla et al. (2008) obtain log [Si]/[Mg] = -0.06 ± 0.05, within 
the F and G dwarf distribution. The GD 40 star+planetary system was probably not 
born with the distinctive Si/Mg ratio observed in the accreted material. One possibility 
is that the GD 40 parent body experienced a substantial differentiation, such that silicon 
was concentrated in a crust, and Mg in a mantle, and that the crust has been subsequently 
lost. 



